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ABSTRACT 

We have obtained deep photometry in two 1° x 1° fields covering the close pair of dwarf spheroidal 
galaxies Leo IV and Leo V and part of the area in between. From the distribution of likely red giant 
branch and horizontal branch stars in the data set, we find that both Leo IV and Leo V are significantly 
larger than indicated by previous measurements based on shallower data. With a half-light radius of 
r/ l =4.'6±0.'8 (206±36 pc) and r/j=2f6±0f6 (133±31 pc), respectively, both systems are now well within 
the physical size bracket of typical dwarf spheroidal Milky Way satellites. Both are also found to be 
significantly elongated with an ellipticity of e ~ 0.5, a characteristic shared by many of the fainter 
(My > —8) Milky Way dwarf spheroidals. The large spatial extent of our survey allows us to search 
for extra-tidal features in the area between the two dwarf galaxies with unprecedented sensitivity. The 
spatial distribution of candidate red giant branch and horizontal branch stars in this region is found 
to be non-uniform at the ~ 3er level. Interestingly, this substructure is aligned along the direction 
connecting the two systems, indicative of a possible 'bridge' of extra-tidal material. Fitting the stellar 
distribution with a linear Gaussian model yields a significance of 4a for this overdensity, a most likely 
FWHM of ~16 arcmin and a central surface brightness of ~32 mag arcsec -2 . We investigate different 
scenarios to explain the close proximity of Leo IV and Leo V and the possible tidal bridge between 
them. Orbit calculations demonstrate that the two systems cannot share the exact same orbit, while 
a compromise orbit does not approach the Galactic center more than —160 kpc, rendering it unlikely 
that they are remnants of a single disrupted progenitor. A comparison with cosmological simulations 
shows that a chance collision between unrelated subhalos is negligibly small. Given their relative 
distance and velocity, Leo IV and Leo V could be a bound 'tumbling pair', if their combined mass 
exceeds 8±4 x 10 9 M Q . The scenario of an internally interacting pair that fell into the Milky Way 
together appears to be the most viable explanation for this close celestial companionship. 

Subject headings: galaxies: individual (Leo IV dSph, Leo V dSph) - Local Group 



1. INTRODUCTION 

In the past four years, the Sloan Digital Sky Survey 
(SDSS) has revealed satellite galaxies of the Milky Way 
(MW) that are up to 100 times fainter than those known 
before, more than doubling the number of known satel- 
lites. Most of these new discoveries are dwarf spheroidal 
galaxies (dSph) that are dark-matter dom inated objects 
(|Martin et~al1l2007t ISimon fc Gehall2007D , yet have sur- 
prising ly complex star formation h istories (Ide Jong et al.l 
2008al) and/or morphologies (e.g. IColeman et al.l l2007a: 
Martin et all l2008at Ide Jong et all l2008b| ) . Thanks to 
these new objects and the uniform sky coverage pro- 
vided by SDSS, it is now possible to probe the faint 
end of the satellite galaxy luminosity function (e.g. 
iKoposov et alJ l2008h . This extended data set has en- 
abled a detailed comparison with galaxy formation mod- 
els, leading to better agreement between models and 
observation, thereby largely solving the 'missing satel- 
lite problem', one of the stumbling blo cks for current 
cosmological models on s mall scales (e.g. IKoposov et al.1 
120091: iMaccio et aDl2009f ). Recent studies of the struc- 
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tural properties of these newly discovered dwarf galaxies 
based on SDSS data have revealed that they 
class, more elongate d than their brighter counterparts 
(|Martin et al.ll2008b[) . indicating that some systems may 
be tidally disrupted by the MW. It is as yet unclear 
whether these objects represent a lower mass limit to the 
galaxy formation process, or whether their stellar bodies 
are remnants of originally larger systems. 

I n this context, the rec ently discove red systems Leo 
IV (IBelokurov et al.l [2007^ and Leo V (|Belokurov et al.1 
2008) are particularly interesting, as they are separated 
by less than 3.0° on the sky with fairly si milar distances 
(154 ± 5 kpc and ~ 180 kpc , respectively; iMoretti et al.1 
120091: IBelokurov et al.|[2008h and radial velocities (10.1 ± 
1.4kms _1 and 60.8 ± 3.1 kms -1 , respec tively, with re- 
spect to the Galactic Stand ard of Rest; ISimon fc Gehal 
I2007t IBelokurov et all 120081), implying that they might 
be related to a sing l e disr upting or disrupted progeni- 
tor. IBelokurov et al.l (|2008h also found the distribution of 
blue horizontal branch (BHB) stars in Leo V to be elon- 
gated and even detected two candidate red giant branch 
(RGB) stars at ~13' North of the system, well beyond the 
estimated tidal radius. F uthermore, with its es timated 
half-light radius of -40 pc (|Belokurov et al.ll2008h . Leo V 
appears to fall in the 'size gap' between gl obular clusters 
and c lassical dwarf spheroidal galaxies (cf. IGilmore et al.1 
12001 . 

We have obtained deep photometry in two 1° x 1° 
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Fig. 1. — Layout of observed fields and color-magnitude diagrams. Left: Black solid lines outline the two 1° X 1° fields that were 
obtained, with the gray lines showing the individual chips. The dashed square outlines the third field that was lost to bad weather. The 
locations of Leo IV and Leo V are indicated. Middle: Color-magnitude diagram of Leo IV within 3' from the center. Overplotted is a 14 
Gyr isochrone with [Fe/H]=— 2.3 and m — Af=20.94. Two selection boxes to select HB and RGB stars are shown with solid outlines, while 
a third selection box used for the maximum-likelihood fitting procedure is indicated with a dashed outline. R ight: As the middle panel, 
but for Leo V. Here the isochrone has m — Af=21.23. In both cases the isochrones are from lDotter et al.l IpOOl l. 



fields, covering both Leo IV and V as well as much of 
the area in between them. These new data enable a 
re-determination of the structural parameters for both 
galaxies, as well as a sensitive search for possible extra- 
tidal material connecting the two systems. The remain- 
der of this paper is structured as follows. In $3] the data 
and data reduction process are detailed. Revised struc- 
tural parameters for Leo IV and Leo V are presented in 
The different ways in which we search for evidence of 
extra-tidal stars in between the two systems is described 
in 21 and in SjS] several scenarios to explain the relation 
between them are discussed. A short summary and our 
conclusions are presented in 

2. DATA 

Deep, wide-field photometry in the Johnson B and 
V filters was obtained for two 1° x 1° fields, using the 
LAIC A imager on the 3.5m telescope at the Calar Alto 
observatory. Four 4kx4k CCDs with a pixel size of 0'.'22 
are positioned such that a mosaic of 4 pointings provides 
a 1° x 1° field of view. Observations were carried out 
during the nights of February 24th and 25th 2009, under 
photometric conditions with seeing ranging from l'/2 to 
l'/8. Total exposure times in B and V were 1950s and 
2600s, divided over single exposures of 650s. The loca- 
tions of the fields are shown in the left panel of Figure [TJ 
a planned central field could not be obtained due to bad 
weather. 

After bias subtraction and flatfielding, the data from 
the individual chips were as trometrically corrected and 
stacked using the SCAMP (jBertinl [2006) and SWARP 
software 4 , developed by Terapix. A chip-by-chip analy- 
sis of the data was necessary because of seeing differences 
between different pointings and PSF distortions. Object 
detec tion was done using SExtractor ([Bertin fc Arnouti 
119961) and subsequent PSF-fitting photometry with the 
daophot package in IRAF. By cross-m atching objects 
with SDSS DR7 (|Abazaiian et alj|2009h and converting 
SDSS magnitudes to B and V using the transformations 

4 SCAMP an d SWARP can be obtained from 
http://tcrapix.iap.fr 



from I Jester et al.l (|2005l ). the photometry was calibrated 
to SDSS. The photometry of all stars is made available 
in Table [TJ which contains the coordinates, magnitudes 
and their uncertainties, and the extinction at the posi- 
ti on of each star interp olated from the extinction maps 
of ISchlegel et al.l JT998). Using these extinction values, 
all stars were individually corrected for the limited fore- 
ground extinction towards these fields {E{B — V) ~ 0.1). 
The final calibrated and dereddened magnitudes will 
hereafter be referred to as Bo and V®. Color- magnitude 
diagrams (CMD) of Leo IV and V within radii of 3' are 
presented in Figure [TJ Artificial star tests to assess the 
completeness in all CCDs were done by placing 4 500 
sources spread over a magnitude range of 17 to 25 in 
each CCD and for each filter and analysing them fol- 
lowing the exact same procedure. The completenesses 
derived are shown in Figure [2] for the Leo IV field and in 
Figure [31 for the Leo V field. Seeing variations cause the 
completeness to be non-uniform, but ~80% completeness 
is reached to at least Bq = Vq ~23.5 over the whole field, 
or approximately 1.5 magnitudes fainter than the SDSS. 

3. REVISED PROPERTIES OF LEO IV AND LEO V 

An accurate distance to Leo IV has been determined 
bv lMoretti et al.l (|2009f) using RR Lyrae stars. Although 
the absolute calibration of our photometry would intro- 
duce uncertainties in the determination of the distance 
to Leo V, the luminosity of its HB stars can be com- 
pared to the luminosity of the HB in Leo IV, providing 
a distance offset with respect to the latter. The stars on 
the horizontal part of the horizontal branch (HB) of Leo 
IV have a mean Bq =21.78±0.06, while for Leo V we 
find Bo =22.06±0.04. Assuming a distance modulus of 
20.94±0.07 (154±5 kpc. iMoretti eTall [20091) for Leo IV, 
this gives a distance modulus of 21.22±0.10, or a distance 
of 175±9 kpc, for Leo V. 

Based on our new photometry, we re-determine the 
structural parameters of the two dwarf galaxies using 
the alg orithm developed and detailed in iMartin et al.l 
(2008bJ). We refer the reader to this paper for details 
but recall that this maximum likelihood (ML) algorithm 
uses the spatial positions of stars that are plausible 
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TABLE 1 

Photometric data 



Chip ID* 


a 


S 


B 




V 




E(B-V) 




(° J2000.0) 


(° J2000.0) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


0.11 


173.33533 


2.05084 


23.302 


0.042 


21.902 


0.014 


0.028 


0.11 


173.23705 


2.05027 


23.603 


0.048 


22.295 


0.021 


0.028 


0.11 


173.28453 


2.05040 


23.531 


0.042 


23.185 


0.047 


0.028 


0.11 


173.29111 


2.04999 


24.160 


0.071 


23.551 


0.074 


0.028 


0.11 


173.27151 


2.05098 


23.059 


0.028 


22.804 


0.032 


0.028 



Note. — The complete table is available online. 
* The first digit of the chip ID indicates the field, where a '0' stands for the Northern 
field containing Leo V and a '1' for the Southern field containing Leo IV; the integers 
following the underscore signify the exposure and the chip number, respectively, both 
running from 1 to 4. 
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Fig. 2. — Completeness as a function of magnitude for all CCDs 
in the 1° xl° field containing Leo IV. The solid line is for the B- 
band, and the dashed line for the V-band. Gray vertical lines indi- 
cate the magnitude limits used for the maximum-likelihood analy- 
sis of the stellar spatial distribution: 23.9 for Bo (solid) and 23.4 
for Vb (dashed). 



20 22 24 20 22 24 20 22 24 20 22 24 
Magnitude 

Fig. 3. — As Fig. [2]but for the field containing Leo V. 




members of a dwarf galaxy (selected by a rough color- 
magnitude cut) to determine the best elliptical, expo- 
nential surface density profile of the satellite, and fit for 
the fore/background stars at the same time. Given the 
positions of Leo IV and V near the edge of the observed 
fields, we use the adapted vers i on of the algorithm pre- 
sented in §3.3 of iMartin et~aT1 (|2009D . This version ac- 
counts for holes and gaps in data sets by iteratively pop- 
ulating un-observed regions with a stellar density that 
follows the best model from the previous iteration. The 
contribution of this stochastic method to the uncertain- 
ties on the model parameters is determined via a Monte 
Carlo approach, in which 100 fits are performed for each 
galaxy. 

The results of the application of the algorithm on the 
LAICA data are listed in Table 1 and show that both 
galaxies are quite elliptical, with axis ratios of 2:1. For 
Leo IV, this value is consistent with the previous mea- 
surement based on SDSS data given its large uncertain- 
ties, bu t the ellipticity-half- light radius covariance (see 
Fig. 9 of lMartin et a l. 2008b) implies a system that is sig- 
nificantly larger than previously measured. As for Leo V, 
previous estimates of its properties relied on a circular 
model and therefore also underestimated its size. With 
the new measurements, both galaxies have half-light ma- 
jor axes of over 100 pc, removing them from the 'size gap' 
that exists between dwarf ga laxies and globular c lusters 
at brighter magnitudes (e.g. iGilmore et afll2007j ). The 
major axes of the systems do not seem in any way cor- 
related with their relative positions. The high elliptic- 
ity for Leo IV is, however, different from the results of 
ISand et all (120091), who find Leo IV to be very round, 
bmce ISand et all 12009) use subgiant and main-sequence 
stars that are below our detection threshold and ignore 
BHB stars, differences could be caused by morphologi- 
cal differentiations between different stellar populations. 
The fact that Leo IV and Leo V lie very close to the 
edge, possibly even partly outside of our surveyed area 
of view might also influence our measurements in some 
unforeseen way. 

Using the number of stars in each dwarf galaxy, which 
is an output of the ML algorithm, and assuming a stel- 
lar population (in this case old, 14 Gyr, and metal-poor, 
[Fe/H]=— 2.3), the total luminosit ies are also de r ived fo l- 
lowing the same procedure as in IMartin et al . (2008b) . 
The resulting luminosities are My = —5.8 ± 0.4 and 
—5.2 ± 0.4 for Leo IV and V, respectively. The revised 
value for Leo IV is 0.8 magnitudes brighter than the value 
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TABLE 2 

Structural parameters of Leo IV and Leo V 



Parameter 


Leo IV 


Leo V 


« (J2000) 


11 32 58.6 ± 1.6 


11 31 08.4 ± 1.6 


S (J2000) 


-00 33 06 ± 54 


±02 12 57 ± 12 


9 (deg) 


-59 ± 9 


-84 ± 13 


e 


0.49 ± 0.11 


0.50 ± 0.15 


rjj (arcmin) 


4 6 +0 - 8 
^■"-0.7 


2.6 ± 0.6 


*h (P c ) 


206^31 


133 ± 31 


D (kpc) 


154 ±5«) 


175 ± 9 


m — M (mag) 


20.94 ± 0.07 °) 


21.22 ± 0.10 


[Fe/H] 


-2.3 ± 0.1 a > 


-2.0 ± 0.2 fc ) 


My 


-5.8 ± 0.4 


-5.2 ± 0.4 


Lv (i Q ) 


1.8 ±0.8 x 10 4 


1.0 ±0.5 x 10 4 


/iy (mag/arcsec 2 ) 


27.5 ±0.7 


27.1 ±0.8 


t)GSR (kms^ 1 ) 


10.1 ± 1.4°) 


60.8 ± 3. l d ) 



Note. ^Moretti et~aTl (f200gh ; ^Walker et~aTI 

(12009T) ; ' HSimon fc Gehal (2007ft : d teelokurov et al.l ||20081 '|; 
all other values based on this publication. 



from lMartin et all ()2008bf ). but they agree to within la. 
It is also statisticall y equivalent to the value measured 
by I Sand et al.l (|2009[ ) from their deeper data. The differ- 
ence is due to the fact that the deeper data reveal that 
Leo IV is actually more elongated and has a larger half- 
light radius than that measured from SDSS data. Leo 
V is found to be alm ost a magnitude bright er than the 
lower limit given bv iBelokurov et al.l (|2008D . who only 
took into account the light inside an aperture of 3'. 

4. A BRIDGE OF STARS BETWEEN LEO IV AND 
LEO V? 

At the Galactic latitudes where Leo IV and V are lo- 
cated, b ~ 50°, the Galactic stellar populations do not 
vary significantly. Hence, in the absence of tidal features 
related to Leo IV and/or Leo V the spatial distribution 
of stars, in any region of the color-magnitude plane, is 
expected to be very close to uniform in the area between 
the two galaxies. We define a coordinate system (X, Y), 
rotated with respect to (a,5), such that Leo IV is lo- 
cated at the origin and Leo V lies along the F-axis. Us- 
ing the color-magnitude selection boxes overplotted on 
the CMDs in Figure [TJ candidate HB and RGB stars at 
the distance of Leo IV and V are selected. Their spatial 
distributions in the (X,Y) coordinate system are shown 
in Figure |4j The magnitude limit of Bq <23.5 ensures 
that no systematic completeness biases affect these dis- 
tributions (see Fig. [3]and[2]). Leo IV and Leo V are 
visible as overdensities at the very bottom and top of 
the field. The flattening of Leo IV is also apparent from 
the contours showing the RGB star density, consistent 
with our structural parameter estimates. However, no 
obvious extra-tidal features are visible. 

As a sensitive test for a possible overdensity of extra- 
tidal stars between the two dwarf galaxies, we divide the 
observed field in 0.2°-wide bins along the A-axis. In the 
following analysis we exclude all data within 5 half-light 
radii of the dwarf galaxies (the ellipses in Fig. 0} . The 
combined density of RGB and HB stars in each bin is 
plotted in panels (a) to (c) of Figure [5] with filled cir- 
cles for the individual mosaics and the complete data 
set. The same is done for MW foreground M dwarf stars 
(selected in the color-magnitude region 21. 0< Bq <23.5, 



1.1< Bq — Vq <1.6) and plotted with open squares. Sim- 
ilar 0.2°-wide bins are defined along the <5-axis and their 
stellar densities are shown in panel (d) . The distribution 
of MW foreground stars is consistent with being uniform, 
both in X and in 8. In contrast, the RGB and HB stars 
are uniformly distributed in S, but not in X, where in 
both fields substructure is visible that resembles a bump 
or overdensity at X ~ —0.1°. Based on basic Poisson 
statistics the significance of this bump is approximately 
3.5<7. Kolmogorov-Smirnov (KS) tests (|Press et al.lll992|) 
indicate that the probability of the distribution of RGB 
and HB stars along the A-axis being drawn from a uni- 
form distribution is only 0.4%, thus rejecting this null 
hypothesis at the ~ 3cr level. Along the 5-axis this prob- 
ability is 56%, and for the MW foreground stars the prob- 
abilities are 87% and 62% respectively, thus all consistent 
with a uniform distribution. Both tests therefore indicate 
the presence of substructure only in the distribution of 
RGB and HB stars and only along the direction parallel 
to a line connecting Leo IV and V. In addition, the KS 
test also shows that the distributions of RGB and HB 
stars along the A-axis in the individual fields are consis- 
tent with being drawn from the same distribution, with 
a probability of 59%. The substructure is therefore not 
isolated to a single field. 

To estimate the significance of this overdensity rigor- 
ously, we performed a ML fit similar to the one used to 
determine the structural parameters of the galaxies. The 
density at any point (X, Y) of the field of view, £(X, Y), 
is here modeled by a flat background, with an ad- 
ditional linear enhancement that is constant along the 
Leo IV/V axis (the F-axis of Fig. [4j and of Gaussian 
form G{X\ A, Xq, a) along the A-axis: 

£(A, Y\A, X , a) = G(X\A, X , a) + E b 

A ( 1 [X- X Q 

= _exp -- 



(1) 



27TCT 



A denotes the amplitude of the Gaussian, Xq its center 
along the A-axis and a its standard deviation, also along 
the A-axis. The probability of the model parameters 
(A, Xq, a), given the location (Aj,Yi) of all considered 
stars in the data set, is then 

p(A, Xo,a\Xi,Yi,V stars i) oc 

C(X tl Y 2i y stars i\A, X , a)p(A, X , a), (2) 

where £(A^, Yi, V stars i\A, Xq, a) is the likelihood of the 
model given the data points and p(A, A , a) is the prior 
probability distribution for the model's parameters. The 
likelihood itself is the product, over all stars, of the model 
defined in equation ([T]), taken on the stellar positions. 
Henceforth 

C{Xi,Yi,V stars i\A,X ,a) = ]J X(A;, Y t \A, X , a) 



n 



.4 



27TCT 



exp 




I (-3) 



The background density, can also be expressed as a 
function of the model parameters by requiring that N*, 
the total number of stars in the data set, corresponds to 
the integral of the density, E, over the field of view, of 
area A: 



N* 



G(X\A,X ,a) dXdY)A' 



(4) 
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FlG. 4. — Spatial distribution of HB (blue dots) and RGB 
(black dots) stars around Leo IV and V, selected using the color- 
magnitude boxes drawn in Figure [T] The coordinates are rotated 
so that Leo IV and V are both at X = and Leo IV lies at the 
origin. Red contours show the RGB star density, and correspond 
to densities exceeding 1.5ct, 3<t, 4<t, and 5ct over the background 
density. The gray ellipses indicate the regions that are neglected 
when looking for extra-tidal material, and the dotted vertical lines 
outline the bins used for the histograms in panels (a) to (c) of Fig. 
1 

The ML fits benefit from a larger sample of stars than 
provided by the selection boxes in Figure ED to constrain 
the fit parameters optimally. Here we therefore use a 
different selection box (the dashed box in Figure [T]) con- 
taining 5 891 stars in the surveyed area but more than 
5rh from the center either of the two dwarf galaxies. The 
prior probability distributions for the model's parameters 
were chosen to be uniform over ranges that are physically 
motivated by a stream of stars originating from the dwarf 
galaxies: the probability on a is uniformly distributed be- 
tween 3' (the typical size of the dwarf galaxies) and 12', 
whereas Xq and A are uniformly chosen between ±0.5° 
(the width of the data set), and ±450, respectively. From 
these, we find that the most probable model has the fol- 
lowing properties: X a = -0.13 ±0.03°, a = 0.11 ±0.02°, 
and A = 172+4Q stars/deg 2 . This corresponds to a back- 



ground density of Ej, = 3224stars/deg 2 for the chosen 
selection box and a total of 275 stars in the overden- 
sity for the selection box and footprint used. Comparing 
the peak amplitude of the overdensity to the number of 
RGB and BHB stars in Leo IV yields a central surface 
brightness of \iy ~32 mag arcsec -2 . The 2-D probability 
distribution functions are shown in Figure [5] for all pa- 
rameters. The significance of the peak in the probability 
distribution function is just over 4<r in all cases, leading 
us to conclude that the data show evidence of a 4cr over- 
density of stars that coincides with the bump visible in 
Figure 

5. DISCUSSION 

All tests — the histograms in Figure 03 the KS tests 
and the ML fits — indicate that there is structure in 
the distribution of stars consistent with being Leo IV or 
Leo V members in the area between the systems. The 
strongest case is provided by the last test, which indicates 
a significance of the overdensity of >4er. Aside from the 
preferred model, there is another local maximum visible 
in Figure |6] that corresponds to an underdensity in the 
data set. This coincides with the drop in density visible 
at X ~ +0.2° in panels (b) and (c) of Figure [SJ The 
limited extent of our survey perpendicular to the possi- 
ble bridge of stars implies that an overdensity can also be 
fit as a region with high background, with an adjacent 
underdensity. However, marginalizing over all parame- 
ters but A (top panel of Fig. [6]) makes it clear that the 
probability of this 'underdensity' is only 2%, whereas the 
probability of an overdensity corresponds to 98%. This 
shows that the structure in the field is more likely caused 
by an overdensity than by an underdensity. This infer- 
ence is based purely on statistics, not taking into account 
that it is not clear what could cause an underdensity of 
stars in the field. 

In their recent analysis, ISand et all (|2009| ) find no ev- 
idence for extra-tidal features emanating from Leo IV. 
Due to their much smaller field-of-view their analysis 
is, however, considerably less sensitive to the presence 
of a possible tidal stream, despite deeper photometry. 
In the region between Leo IV and Leo V, where any 
putative bridge would be found, our surve y covers an 
area s even times larger than that available to ISand et "all 
(2009). At a given depth, this leads to a seven- fold in- 
crease in the number of stellar tracers, in our case RGB 
and B HB stars, th a t can be found. The higher sensitivity 
of the ISand et ail f|2009T ) data only partly compensates 
for this, and only adds fainter stars in a region of the 
CMD where cont amination by faint galaxies is high. The 
detection limit of ISand et al.l (|2009l ) for tidal streams of 
H g <29.8 mag arcsec -2 is insufficient to detect the stellar 
bridge, which has an estimated central surface brightness 
of fiy ^32 mag arcsec -2 . It is, however, interesting to 
note that one of t he significant 'nuggets' of stars seen by 
ISand et "all (|2009[ ) to the north-west of Leo IV lies exactly 
on top of the location of the stellar bridge, as indicated 
by our ML fits. 

The spatial proximity of Leo IV and Leo V and the in- 
triguing possibility of extra-tidal material between them 
elicits the question of whether these two systems are 
somehow related and/or interacting with each other. 
Here we consider the following different scenarios: 
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Fig. 5. — Stellar densities in 0.2°-wide strips. In panels (a), (fe) and (e), the strips are at constant X , parallel to the straight line 
connecting Leo IV and V. From left to right, the panels correspond to the field containing Leo IV, the field containing Leo V, and the 
sum of both fields. Elliptical regions corresponding to 5 half-light radii around Leo IV and Leo V are avoided (Fig. |4j. In all panels, filled 
circles are for stars selected with the RGB and HB color-magnitude selection boxes (Fig. [1} and open squares for foreground M dwarf stars. 
Vertical error bars show the Poisson errors and the horizontal error bars the width of each strip. In panel (d) the strips are at constant 
declination 5, almost perpendicular to the straight line connecting Leo IV and V. 



• Leo IV and Leo V might be condensations of stars 
in a much larger stellar stream, even though the 
density contrast between the condensations and the 
underlying stream would be very high. The appar- 
ent width of the extra-tidal substructure in the can- 
didate RGB and HB stars of — 15 arcmin is similar 
to the sizes of Leo IV and V, and is therefore con- 
sistent with this scenario. In this case they would 
be expected to closely follow the same orbit, which 
should come close enough to the Galactic center to 
cause the tidal disruption of a common progenitor. 

• Leo IV and Leo V might have fallen into the 
Milky Way together as members of a group of 
dwarf galaxies, a process exp ected to be com- 
mon i n a ACDM cosmo logy (e.g. lD'Onghia fc Lake! 
[20081 : iLi fc Helm! [20081. or as satellites of a larger 
halo, as has been su ggested for Segue 1 and 2 
(Bclokuro v et al.ll2009h . Regarding the latter case, 
it should be noted that the revised half-light radii 
of Leo IV (-200 pc) and Leo V (-130 pc) make 
them much larger systems than Segue 1 and 2 and, 
in terms of size, clearly more similar to other dSphs. 
If they do constitute a pair with a common origin, 
they would still follow the same orbit around the 
Milky Way, but with a possibly large velocity dis- 
persion, as they orbit a common center of mass. 
In this scenario, extra-tidal stars might be the re- 
sult of interactions between the two galaxies, rather 
than with the Milky Way. 

• A final option would be that Leo IV and Leo V were 
originally unrelated satellites of the Milky Way, but 
had a chance encounter. A near head-on collision 
might be able to cause an interaction strong enough 
to disrupt their stellar bodies. 

5.1. A common orbit 

If Leo IV and Leo V are tied together and follow the 
same orbit around the Milky Way, their common orbit 
should be easily determinable from simple physical ar- 
guments. In the following, we assume that the Galac- 
tocentric properties of the galaxies are interchangeable 
with their heliocentric properties, a reasonable assump- 
tion given the large distances to the systems compared to 



the solar radius around the Milky Way. If the two dwarf 
galaxies lie along the same orbit, they should actually 
share the same orbital energy and angular momentum. 
The latter requirement can first be used to link together 
the tangential velocities of the two galaxies in the plane 
of the sky, vt.i and 1^5, and their distances, r± and r$, 
respectively for Leo IV and Leo V: 1^5 = Vt,4 r 4/ r 5- We 
then invoke energy conservation along the Leo IV/ V or- 
bit to obtain 



(«tot,4) - (Vcsca) = («tot,5) - (' 



x2 
■'esc, 5/ 



(5) 



where Vtat,i — y i v t.i) 2 + {vr.i) 2 is the total velocity of 
Leo i, v r ,i its radial velocity with respect to the Galactic 
Standard of Rest, and w CS c,i the escape velocity in the as- 
sumed model of the Milky Way potential, at the location 
of Leo i. For Leo IV, the tangential motion of the orbit 
that links Leo IV and Leo V is therefore defined by: 

via = 



\ 



~ (6) 



In order to determine 1^4 and 1^5 with their associ- 
ated uncertainties, we follow a Monte Carlo scheme and 
randomly generate the distances and radial velocities for 
each of the two dwarf galaxies from Gaussian distribu- 
tions defined by their observed properties, as listed in Ta- 
ble [2j We then determine the escape velocities required 
by equation (EH) for the Milky Way model adopted by 
iPaczv hski (1990), after removing unphysical cases where 
one of the two galaxies is not bound to the Milky Way. 
The median and central 68.3% of the resulting distribu- 
te km s -1 and 

1 



tions after 10000 trials yield v t ,4 = 239_ 6 



Vt,5 



211±fg kms" 



These large values, combined with the small radial ve- 
locities of the dwarf galaxies, make it impossible to find 
a viable orbit that reproduces the observed properties 
of both systems at the same time. None of the orbits 
drawn in Figure [7] goes through both points in all pan- 
els, demonstrating that, whether Leo IV or Leo V is used 
as a starting point, neither the velocity nor the distance 
along the orbit has a strong enough gradient to replicate 
the observed properties of the other dwarf galaxy. This 
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Fig. 6. — Probability distribution functions from maximum-likelihood fits of a linear Gaussian to the stellar distribution between Leo 
IV and Leo V. Top panel: probability distribution function of the model parameter A marginalized over the two other parameters. The 
region of negative A, shown in grey, yields a small probability of only 2%, compared to 98% for the region of positive A; the probability 
of the best model (indicated with a filled circle in all panels) is a factor ~10 4 higher than that of a uniform distribution (A=0), which 
corresponds to a ~4cr significance. Bottom panels: for all combinations of the three fit parameters (A, Xo and cr) the two-dimensional 
probability function is shown, marginalized over the third parameter. Here, contours represent drops of 50%, 90%, 99% and 99.9% of the 
probability with respect to the best model. 



analysis thus shows that Leo IV and V cannot be follow- 
ing the exact same orbit. We can calculate a compromise 
solution, with a starting point for the orbit integration 
that has the mean properties of the two galaxies (the cen- 
tral line in Figure [7]) . This orbit does not reproduce the 
velocities and distances of the two systems and remains 
at large Galactocentric distances, with an apocenter of 
244 kpc and a pericenter of 158 kpc. We conclude that 
Leo IV and Leo V are very unlikely to be condensations 
in a larger stellar stream, for two reasons: no viable single 
orbit exists, and the compromise solution does not come 
close enough to the Galactic center for a single progenitor 
to have been tidally affected. 

5.2. A 'tumbling pair' 

Rather than being remnants of a progenitor disrupted 
by the gravitational tides of the Milky Way, the two 
dwarf galaxies could be a gravitationally bound pair in- 
teracting with each other. This means that they would 
follow a 'common' orbit, but with a large velocity disper- 
sion, as they orbit their common center of mass. Such a 
'tumbling pair' of dwarf galaxies would be similar to the 



Large and Small Magellanic Clouds, albeit much more 
scaled-down. For the two systems to be gravitationally 
bound, the total kinematic energy of the system must 
be less than th e total gravitation al potential energy. As 
pointed out bv lDavis et al.l (|1995l ). the two-body Newto- 
nian binding criterion can be written as 



> 



2G sin 2 a ' 



(?) 



where M sys is the total mass of the system, 7? the dis- 
tance between the two objects, vi os the line-of-sight rela- 
tive velocity, and a the angle between the axis of the two- 
body system and the sky. Inserting the values for these 
parameters (u; os = 50.7 ± 3.4 km s _1 , R — 22 ± 10 kpc, 
a = 70° jlf 8 ; see Tableland references therein) yields a 
lower limit for the total mass of M sys > 8 ± 4 x 10 9 Mq, 
or half of that for each individual dwarf. 

Masses inferred from the stellar velocity dispersions in 
faint dSph galaxies are typically much lower, but only 
probe the inner few hundred parsec of their dark mat- 
ter halos. Based on dynamical modeling of dSphs em- 
bedded in cosmologically motivated dark matter halos, 
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Fig. 7. — Orbits for Leo IV and V, derived assuming common 
energy and angular momentum for the two dwarf galaxies. The 
dotted, dashed and solid lines correspond to orbits with initial con- 
ditions determined using Leo IV, Leo V or their mean properties, 
respectively. The panels show the evolution of Galactic latitude 
(top), heliocentric distance (middle) and radial velocity with re- 
spect to the Galactic Standard of Rest (bottom) as functions of 
Galactic longitude. Properties of Leo IV and V are indicated by 
filled and open circles, respectively. 




100 
D .ub-.ub [ k P c ] 



1000 



Fig. 8. — Cumulative probability distribution for the distance 
between two dark matter subhalos (with a current mass larger than 
5 X 10 7 M Q ) in the last 3.4 Gyr (see text for details). 



iPenarrubia et al.l ()2008[ ) derive masses of up to several 
times 10 9 M Q for Local Group dSphs. The masses re- 
quired for Leo IV and Leo V to be a bound pair are 
therefore possible, but imply high mass-to-light ratios of 
M/Ly=10 5 in solar units. 

5.3. A close encounter 

A third option would be that Leo IV and V are un- 
related systems that randomly collided in the Milky 
Way halo, thus provoking an interaction. Given that 



the objects have a stellar extent of ~500 pc, we esti- 
mate that they would need to come to within a few 
kiloparsecs of each other with a relative velocity of less 
than ~ 100 km s -1 for tides to affect their stellar com- 
ponents. We used results from N-body simulations in 
order to compute the probability of such an encounter. 
We started from the four Milky Way-like dar k matter 
(ACDM) halos presented in lMaccio etaLI (|2009t ) and, for 
each subhalo within 250 kpc of the host galaxy's center 
and with a bound mass today M su b > 5 x 1O 7 M , we 
determined the distance to any other subhalo, D su b~ S ub, 
in the last 3.4 Gyr. Figure shows the cumulative prob- 
ability distribution for D su b- su b in the four Milky Way 
analogue simulations. The probability of two halos com- 
ing to within 5 kpc of one another is as low as a few 
times 10 - °. In addition, this number can only be an 
upper limit for a galaxy-galaxy interaction, since only 
« 30% of dark matter halos with M ?u b > 5 x 10 7 M Q are 
expected to host stars (jMaccio et al.ll2009l ). We conclude 
that a collision between Leo IV and Leo V is not a likely 
explanation for a possible interaction. 

6. CONCLUSIONS 

Given that Leo IV and Leo V form a close pair on the 
sky, in distance, and in radial velocity, it has been argued 
tha t they might be interac ting or even share the same or- 
bit (jBelokurov et al.lfeOOSl ). We have obtained new, deep 
photometry in two 1° x 1° fields containing both dwarf 
spheroidal galaxies, re-derived their structural proper- 
ties, and searched for extra-tidal stars. Both Leo IV and 
Leo V are considerably larger than previously thought, 
with sizes of 206 pc and 133 pc. These sizes are simi- 
lar to those of other recently discovered dSphs such as 
Bootes I or Ursa Ma jor II (respectively 242 and 140 pc, 
iMartin et al.1 l2008bT ) and t he smaller 'classical' dSphs, 
for examp le Draco (221 pc. [Martin et al.ll2008bf ) or Leo 
II (185 pc. lColeman et al]l2007b| h~ Both are also strongly 
elongated, with axis ratios of 2:1, which i s the case for 
many of the fainter [My > —8) dSphs (|Martin et al.l 
2008bj). We note, however, that the proximity of both 
dwarf galaxies near the edge of our survey might have 
some influence on these measurements. 

The combination of deep photometry and a large area 
coverage of our survey provides an unprecedented sen- 
sitivity to low surface brightness features in the region 
between Leo IV and Leo V. Analysis of the spatial dis- 
tribution of candidate RGB and HB stars using den- 
sity histograms and Kolmogorov-Smirnov tests reveals 
the presence of a significant substructure in the observed 
fields. This substructure is detected at the ~ 3er level, 
but only along the direction connecting the two systems. 
Maximum-likelihood modeling of the stellar distribution 
with a linear Gaussian confirms that the spatial struc- 
ture can be fit with an overdensity connecting the two 
systems. With this method, the significance of the over- 
density is just over 4c, and has a peak surface brightness 
of nv — 32 mag arcsec~ 2 . 

The reason for the close proximity of Leo IV and Leo V 
is still an enigma, and we have considered a few possible 
explanations. Orbits derived for the objects using their 
positions and velocities and by requiring the pair to share 
common values for energy and angular momentum show 
that it is not possible for both to be on the exact same or- 
bit under reasonable assumptions. Furthermore, a com- 
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promise orbit would not approach the Galactic center 
closer than ^160 kpc, as also noted bv iBelokurov et all 
( 2008). Hence, it is highly unlikely that Leo IV and Leo 
V are condensations in a low surface brightness stellar 
stream resulting from the tidal disruption of a common 
progenitor. The fact that the elongations of the two sys- 
tems are not aligned might be a further indication of this. 
A large velocity dispersion around a mutual orbit, for ex- 
ample in the case of a 'tumbling pair' of dwarf galaxies, 
is only viable if Leo IV and Leo V have masses of at 
least a few times 10 9 M©. This would imply extreme 
y-band mass-to-light ratios of >10 5 in solar units, but 
recent dynamical simula tions of dwarf spheroidal galaxies 
(|Peharrubia et all I2008D show that the required masses 
are possible. If Leo IV and Leo V are indeed a bound and 
internally interacting pair, this could account naturally 
for their unaligned elongations and the offset of the ap- 
parent 'stellar bridge' from the straight line connecting 



them. Finally, analysis of Milky Way-like ACDM halos 
reveals that the probability of two satellites colliding is 
negligibly small. This is consistent with the findings of 
IBelokurov et al.l (|2008[ ) who calculate a <1% probability 
for such a close association happening by chance. There- 
fore, the scenario in which Leo IV and Leo V are a bound 
pair of dwarf galaxies, orbiting and interacting with each 
other, appears to be the most viable explanation for this 
close celestial companionship. 
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